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Abstract. - The detection of clusters of galaxies in large surveys plays an important part in extra- 
galactic astronomy, and particularly in cosmology, since cluster counts can give strong constraints 
on cosmological parameters. X-ray imaging is in particular a reliable means to discover new clus- 

Oters, and large X-ray surveys are now available. Considering XMM-Newton data for a sample 
of 40 Abell clusters, we show that their analysis with a Kolmogorov distribution can provide a 
distinctive signature for galaxy clusters. The Kolmogorov method is sensitive to the correlations 
' in the cluster X-ray properties and can therefore be used for their identification, thus allowing to 

| search reliably for clusters in a simple way. 

o 

C3 



> 

in 

On 



O 



X 



Introduction. — The reliable determination of galaxy 
clusters as entities is crucial not only for the studies 
of their individual structure and evolution, but in the 
broader context of the formation of large scale struc- 
tures in the universe. Various methods to define the clus- 
ter membership of galaxies have been developed, usually 
based on limited information such as line-of-sight veloci- 
ties or photometric redshifts (ID), and coordinates (2D). 
X-ray imaging has provided an important means to de- 
tect clusters as diffuse X-ray sources and to define cluster 
properties by allowing to derive their overall gravitational 
potential, including dark matter, under the assumption 
that the X-ray emitting gas is in hydrostatic equilibrium 
in the cluster gravitational potential well dominated by 
dark matter. 

In the present paper we apply the Kolmogorov stochas- 
ticity parameter method and distribution [TJ[2] to analyze 
X-ray data, in order to reveal possible signatures which 
can allow to identify a galaxy cluster. This approach quan- 
tifies the degree of randomness for given sequences of num- 
bers and has been applied to various systems appearing in 
number theory [3JH] . We have also applied this method to 
the cosmic microwave background (CMB) radiation sig- 
nal [5j[8]. It allowed to identify non-Gaussianities in the 



( a ) Based on observations obtained with XMM-Newton, an ESA 
science mission with instruments and contributions directly funded 
by ESA Member States and NASA 



CMB such as the Cold Spot in the maps of the Wilkin- 
son Microwave Anisotropy Probe (WMAP). The behavior 
of the stochasticity parameter supports the void nature 
of that anomaly, and another non-Gaussianity, the North 
Cold spot, has also been identified in the north sky by 
the same method [6]. This approach was also efficient 
at detecting point sources in the WMAP maps, including 
the prediction [7] of gamma-ray sources (quasars, blazars) 
which have later been identified by the Fermi satellite [9]. 

Much has been learnt these past ten years on 
the properties of X-ray clusters, based on the 
data acquired by the XMM-Newton, Chandra and 
Suzaku satellites, and the corresponding references 
are too numerous to be quoted here (for overview 
see e.g. NASA and ESA sites for X-ray missions 
http://www.nasa.gov/missions/index.html; 
http://www.esa.int/esaSC). In particular, with 
its large collecting area, XMM-Newton has allowed not 
only to draw accurate emissivity maps, but also to com- 
pute hardness ratio, temperature and metallicity maps of 
the X-ray gas (as first obtained by [10j[TT]). These maps 
give a much deeper insight on the cluster properties, in 
particular on the cluster merging history. They have 
shown unambiguously that even clusters with X-ray 
emissivity maps showing a relatively relaxed appearance 
could in fact be undergoing one or several mergers (see 
e.g. [H] and references therein). Several cluster surveys 
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based on XMM- Newton data are available, the first 
being those of [T3HT5] . The XMM- Newton archive is now 
extremely large, and a simple method applied to detect 
clusters throughout this archive would be a very powerful 
tool to obtain a large and homogeneous sample of clus- 
ters. Note that presently very large and homogeneously 
detected cluster catalogs are also being obtained, based 
on the Sunyaev-Zel'dovich effect, both with ground based 
facilities such as the South Pole Telescope [16] or the 
Planck satellite [T7] . 

The aim of the present study is not to analyse the prop- 
erties of individual clusters, but to show that the Kol- 
mogorov descriptor, considered as a quantifier of corre- 
lations in datasets, can be used for the identification of 
galaxy clusters based on X-ray data. For this particular 
aim we will use a sample of Abell clusters with available 
XMM-Newton data. 

The paper is structured as follows: we briefly describe 
our method in Section 2, we present our sample of X-ray 
clusters and analysis together with our results in Section 3. 
Conclusions are drawn in Section 4. 

The method. — The definition of the Kolmogorov 
stochasticity parameter [TJ[2] is given for a finite random 
sequence of real numbers x±,X2, . . . ,x n so that the values 
of x n are sorted in an increasing manner: x\ < x^ < • • • < 
x n . 

One can define an empirical distribution function: 

C n (X) = (number of the elements Xi < X), 

where 

( 0, X < x x 
C n (X) = l k/n, x k <X<x k +! (1) 
[ 1, x n <X 

The theoretical distribution function C(X) is then: 

C(X) = n • (probability of the event x < X). 

The stochasticity parameter A n for a sequence of n val- 
ues of x is defined as: 

A n = su Px \C(X) - C n (X)\/yfii. (2) 

According to Kolmogorov's theorem, A n is a random 
number with an empirical distribution function 

^> n (A) = (number of elements X n < A), 

uniformly converging to $(A) at n — )• oo for any continu- 
ous distribution C(X) 

*„(A) *(A), (3) 

where 

+oo 

$(A)= (- 1 ) fe e" 2fc2A2 , -oo<fc<+oo. (4) 

k— — oo 



The distribution function 3>(A) varies monotonically 
from $(0) = to $(oo) = 1. 

According to the definition, for large enough n and ran- 
dom sequence x n the Kolmogorov stochasticity parame- 
ter A n will have a distribution close to 3>(A). If the se- 
quence is not random, the distribution will be different. 
So 3>(A) denotes a quantitative measure for the random- 
ness of a sequence. Arnold has shown the informativity of 
this method already for sequences of n=15 [2]. 

Analysis and results. — Our analysis is based on 
a sample of 40 Abell clusters with data available in the 
second XMM-Newton Serendipitous Source Catalogue, 
2XMjVL released on April 15, 2010: the 2XMMi-DR3 cat- 
alogued in the [0.2-12] keV band, of median flux 2.5 10~ 14 
erg cm -2 s _1 . The X-ray images of the clusters represent 
datasets of X-ray fluxes for pixels with given coordinates. 
Table 1 shows the relevant data for the 40 clusters; the 
first 17 have redshifts lower than 0.03, although no noti- 
cable redshift dependence have been noticed at the study 
below. For illustration in Figure 1 we represent the X- 
ray image for a particular cluster, A2870, and in Figure 
2 the cuts of its X-ray flux distribution over the Galactic 
/ and b coordinates. 

We then create Gaussian mocks of the clusters, namely, 
the mock of a given cluster is obtained as a Gaus- 
sian distribution of the flux with the median and 
standard deviation defined by the real data of its 
X-ray pixelized distribution; the cuts along two axes 
are also shown in Figure 2. Since the Gaussian mocks of 
the clusters cannot possess the information on the clus- 
ter potential carried by the real X-ray images, they can 
serve as a calibration for the correlations in the real data, 
i.e. what we have to quantify using Kolmogorov's method. 
Comparing the plots in the left and right columns in Fig- 
ure 1, we see that although the real cluster data and the 
Gaussian ones exhibit certain discrepancies, they do not 
differ radically, and in any case they cannot be distin- 
guished reliably only via such a comparison. 

We then apply the Kolmogorov technique to the X- 
ray data for each cluster. Namely, we estimate the Kol- 
mogorov stochasticity parameter A and the function <£(A) 
for the X-ray images, and do the same for the Gaussian 
images of each cluster. The maximal and minimal values 
of the resulting functions are shown in Fig. 2 for all 40 
clusters, as well as for their Gaussian mock clusters. We 
see that for real X-ray images of clusters the function <1> 
outside the center remains 1. On the other hand, its be- 
havior for the Gaussian images is entirely different, i.e. far 
smaller at all radii. To test the dependence on the number 
of pixels in Fig. 3 we show these distributions for 4 clus- 
ters with large scatter in the pixel numbers: A0150 - 60; 
A3526 - 248; A0576 - 315; A2052 - 324. It is clear that, the 
higher is number of pixels, more outlined is the difference 
between the real clusters and the mock Gaussian ones. 



http:/ /xmmssc- www. star. le.ac.uk/Catalogue/2XMM/ 
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Table 1: Sample definition for the 40 clusters: running num- 
ber, cluster name, redshift, galactic coordinates (of centers), 
number of pixels for which X-ray fluxes are given; the first 17 
clusters have redshifts lower than 0.03, the others have a higher 
redshift. 



N Abell z 1, deg b, deg Pixel 

97 
183 
116 
85 
86 
102 
758 
324 
99 
45 
45 
103 
73 
169 
248 
162 
285 
83 
100 
159 
77 
60 
65 
127 
217 
217 
315 
94 
73 
265 
113 
98 
112 
143 
145 
88 
111 
78 
393 
113 



Conclusions. — We analyzed a sample of 40 Abell 
X-ray clusters with available XMM- Newton data. By ap- 
plying the Kolmogorov method to the X-ray data and to 
Gaussian simulated images, we revealed a clear difference 
in the behavior of Kolmogorov's function <£> for real clus- 
ters and for their Gaussian simulated images. The X-ray 
image of a real galaxy cluster carries correlations reflect- 
ing the gravitational potential of the cluster and hence, 
the physical conditions of the hot gas responsible for the 
X-ray emission. Such correlations are not present in the 




Fig. 1: X-ray image of the cluster Abell 2870 obtained by 
XMM-Newton. The circle has a radius of 0.3°; the '0' of the 
scale bar below the figure corresponds to a flux of 41.3 in units 
of 10 -15 erg cm -2 s _1 , the minimal and maximal values of 
the bar are 7.3 and 405.6, respectively, in the same units. 



randomly redistributed simulated Gaussian clusters. Al- 
though the original X-ray clusters and their simulated 
Gaussian versions do not show radically different behav- 
iors in their flux distributions, a difference, however, is 
clearly outlined quantitatively in the behavior of the Kol- 
mogorov function <£. For the real X-ray images of the 
clusters the function $ tends to 1 outside the cluster cores, 
while for the Gaussian mock clusters it remains smaller at 
all radii. 

Thus, the Kolmogorov method applied to X-ray images 
enables to distinguish a galaxy cluster from a random X- 
ray flux distribution, thus opening new possibilities in the 
detection of galaxy clusters. 
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Fig. 2: Distribution of the pixel X-ray flux along the / (top) 
and b (bottom) axes for A2870 (left) and for the simulated 
Gaussian distribution of the X-ray data parameters (right). 
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Fig. 3: Maximal and minimal values of the Kolmogorov func- 
tion vs the cluster radius (in degrees) for X-ray data for 40 
galaxy clusters, as well as the same for corresponding Gaussian 
mock clusters (in blue). 
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Fig. 4: Kolmogorov function vs the radius for 4 clusters with 
large scatter in the pixel numbers (along with their Gaussian 
versions): A0150 - 60; A3526 - 248; A0576 - 315; A2052 - 324. 



p-4 



